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We use a non-Markovian approach to study the decoherence dynamics of a qubit (with and 
without measurement) in either a low- or high-frequency bath modeling the qubit environment. 
This approach is based on a unitary transformation and does not require the rotating-wave 
approximation. We show that without measurement, for low-frequency noise, the bath shifts the 
qubit energy towards higher energies (blue shift), while the ordinary high-frequency cutoff Ohmic 
bath shifts the qubit energy towards lower energies (red shift). In order to preserve the coherence 
of the qubit, we also investigate the dynamic of qubit with measurement (quantum Zeno effect) 
in this two cases: low- and high-frequency baths. For very frequent projective measurements, 
the low-frequency bath gives rise to the quantum anti-Zeno effect on the qubit. The quantum 
Zeno effect only occurs in the high-frequency cutoff Ohmic bath, after considering counter-rotating 
terms. For a high-frequency environment, the decay rate should be faster (without measurements) 
or slower (with frequent measurements, in the Zeno regime), compared to the low-frequency bath 
case. The experimental implementation of our results here could distinguish the type of bath 
(either a low- or high-frequency one) and protect the coherence of the qubit by modulating the 
dominant frequency of its environment. 
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I. INTRODUCTION 

There is considerable interest in low-frequency noise (see, e.g., the review in Ref. [jy] and references therein), because 
this type of noise limits the coherence of qubits based on superconducting devices such as flux or phase qubits Q, 0| . 
Also, the dephasing of flux qubits is due to low- frequency flux noise with intensity comparable to the one measured 



in dc SQUIDs (e.g., Refs. [l,|J,|5|). Under certain conditions, noise can enhance the coherence of superconduting 
flux qubit 0. There are also several models (e.g., Refs. for the microscopic origin of low-frequency flux noise 

in Josephson circuits. Therefore, the study of low-frequency noise has becomes very important for superconductor 



pbi.. a Moreover, the <,„a,rt,„. Ze„o effect b,en proposed . strategy to protect colrercncc M, e„t.ng,en,e„t 
[ll[E3|i to control thermodynamic evolutions pLsJ. The quantum Zeno effect and anti-Zeno effect have been widely 
discussed [3] ■ So it is an interesting topic to investigate the quantum Zeno effect of a qubit coupled to a low- 
frequency bath. 

The description of low- frequency noise fl| (such as 1 /f noise) is complicated by the presence of long-time correlations 
in the fluctuating environment, which prohibit the use of the Markovian approximation. In addition, the rotating 
wave approximation (RWA) is also unavailable in an environment with multiple modes. In the case of a time-depend 
external field for a qubit coupled to a thermal bath, Kofman and Kurizki developed a theory [igI which considers 
the counter-rotating terms of the fast modulation field through the negative-frequency part G{lu) (lu < 0) in the 
bath-correlation function spectrum. The dynamics in two significant models (the spin-boson model with Ohmic bath, 
and a qubit coupled to a bath of two level fluctuators) have been calculated within a rigorous Born approximation 
and without the Markovian approximation Refs. 17 1 and [3| describe the structure of the solutions in the 



complex plane with branch cuts and poles. 

Here we present an analytical approach, based on a unitary transformation. We use neither the Markovian approx- 
imation nor the RWA in order to discuss the transient dynamics of a qubit coupled to its environment. This method 
has already been used [l^ to study the decoherence of the Ohmic bath, sub-Ohmic bath, and structured bath. In this 
paper, we calculate the coherence dynamics of the qubit (with and without measurement) respectively in two kinds 
of baths. Besides producing an energy shift, the environment can change the decay rate of the qubit. To preserve the 
coherence, we also investigate the decay rate of the qubit subject to the quantum Zeno effect. The low-frequency noise 
uses a Lorentzian-type spectrum with the peak of the spectrum in the low-energy region and for the high-frequency 
noise we choose the ordinary Ohmic bath with Drude cutoff. 

Our results show that for Zow-frequency noise, the qubit energy increases (blue shift) and an anti-Zeno effect takes 
place. For a /ii(;/i-frequency cutoff Ohmic bath, the qubit energy decreases (red shift) and the Zeno effect dominates. 
For a high-frequency environment, the decay rate should be faster (without measurements) or slower (with frequent 
measurements, in the Zeno regime), compared to the low- frequency bath case. The experimental implementation of 
our results here could distinguish the type of bath (either a low- or high-frequency one) and protect the coherence of 
the qubit by modulating the dominant frequency of its environment. The coherence dynamics without measurement 
also indicates that the coherence time of the low-frequency noise is longer than the high-frequency noise, which 
demonstrates the powerful temporal memory of the low-frequency bath and stem from non-Markov process of qubit- 



bath interaction within the bath memory time 20| . 
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II. METHOD BEYOND THE RWA BASED ON A UNITARY TRANSFORMATION 



We describe a qubit coupled to a boson bath, modeling the environment, by the Hamiltonian 



H 



1 



-Act. 
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gk{al + ak)<Tx + ^ ujka^ak- 



(1) 



The boson bath and the spin-i (fluctuators) bath will have the same dissipative effect on the qubit at zero temperature, 
T = 0, if both baths have the same correlation function 21|. However, for finite temperatures, the spin bath has a 



smaller effect on the qubit because of the likely saturation of the populations in the spin bath. Here, for simplicity, 
we only consider the case of zero temperature. That is, our results are also applicable to a fluctuator-bath at zero 
temperature. Our approach is based on a unitary transformation and can be used for different types of environmental 
baths. Below we give detailed derivations for the low-frequency noise case. 
The spectral density J(a;) of the environment considered here is given by 



Ji^) = '^9kSi^^-^^k) = 



2auj 



(2) 



where A is an energy lower than the qubit two-energy spacing A and a describes the coupling strength between the 
qubit and the environment. Choosing A as the energy unit, a/A^ is a dimensionless coupling strength. When w > A, 
J{uj) ^ l/w, corresponding to a 1// noise. 

We apply a canonical transformation to the Hamiltonian H: H' = exp [S] exp [— S*] . Further explanations on 
the validity of the transformation can be found in the Appendix. Note that J2 k 9k^k'^ / i'^'^k) > in Eq. (|A.7p . so 
< exp [— g^£,k'^ / (2w^)] < 1. Then the solution of -q will be in the region from to 1. Actually, the existence and 
uniqueness of the solution of rj in Eq. (IA.7I) can be used as a criterion for the validity of our method. The parameter 
•q can be regarded as a rcnormalization factor of the energy spacing A and is calculated as 



TrAry A 



rj = exp < a- 



r/2A2-h(A2- r/2A2)log 



)) A 



(A2-|-r;2)2 



(3) 



Obviously, 77 is determined self-consistently by the above equation. 
Thus, the effective transformed Hamiltonian can be derived as 



k k 



where 



9k£.k 



(4) 



(5) 



Comparing H' in Eq. (j4]) with the ordinary Hamiltonian in the rotating- wave approximation (RWA): 

^^RWA = + ^ Wfefl+afc + ^ Y(a+cr_ -I- afecr+), (6) 

k k 

one can see that the unitary transformation plays the role of renormalizing two parameters in the Hamiltonian, i.e., 
the energy spacing A is renormalized: 



A^T^A; 



(7) 
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and the coupling strength 5fc/2 between the qubit and the bath is renormahzed: 



9k_ 
2 



gk 



2f] A 



(8) 



Below we study the decoherence dynamics of the qubit and the quantum Zeno effect using the transformed Hamiltonian 
H'. 

A. Non-measurement decoherence dynamics 



We diagonalize the transformed Hamiltonian H' in the ground state \g) = |t) \0k) and lowest excited states, ||) |Ofe) 
and It) |lfc) , as 



h' =-lrjA\g} {g\ + Y,E\E} {E\ , 



(9) 



where |t) and \l) are the eigenstates of a^, i.e., |t) = It) > f^z It) = — It) and \nk) denotes the state with n bath 
excitations for mode k. The state \E) is 



with x{E) = l + Ek 



\E)=x{E) |t)|0,)+^2;fe(i?) |t)|U 

k 

-1/2 



(10) 



and yk{E) = g+^2/2-c., ^(-^)- 
Here we calculate the dynamical quantity {ax{t)) , which is the analog of the population inversion {az(t)) in the 
spin-boson model [22ij. Since the coupling to the environment will be "always present" in essentially all physically 
relevant situations, a natural ground state is given by the dressed state of the two-level qubit and bath. Therefore, 
the ground state of H is exp[— 5] |t) \0k) , under counter-rotating terms, and the corresponding ground-state energy 
is —1] A/2, then the ground state of H becomes |t) |Ofc) with the identical ground-state energy —rj A/2. We now 
prepare the initial state, which is also a dressed state of the qubit and bath, from the ground state as 



|^(0)) = ^i±|^exp[-^] n)\Ok) 



(11) 



Then the initial state in the transformed Hamiltonian becomes \ip{0)} = (|t) + It)) |0a;) /V2, which is the eigenstate 
of ax- Starting from this initial state, we obtain j23 | 



{<jx{t)) = Ti-B{mw^\m) 

= Trs (^-(0) \exp[iHt](Tx exp[~iHt] \ ^(0)) 



E 



= ^Y^xiEfexp 

) 

dE exp 
dE exp 



-^^E + !l^H 



Ij^xiE^exp 



^^E+!1^H 



1 

Ani 



1 



-iE t 



iE t 



E' + iO+ - ujk 



Re 



1 

27ri 



exp 



-iEt 



E - iO+ - ujk 



E'-VA-Ekw:^.. 



-dE 



(12) 
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Here we denote the real and imaginary parts of J^k I ('^ ~ "^fc ± *0"'") as RiiS) and =Fr(a;), respectively. It follows 
that 



"~ — i^k J 

k 

^ A 9 9 log |W 



[uj — uj'){ijj' + ryA) 



2 



(77A + a;)2(A2 +a;2) 
2Q!A^r/' 



2 2 7rA[A2+77A (-??A + 2a;)] - 2 [A2(2r7 A - cj) + ry^A^ a;]log|A| 



2(A2+7y2A2)2(A2+ Cj2) 

2 A + c^) (A2 + A2r?2) ^ [27?3A3 + (-A2 + r?2A2) u] log ^ A| 

+ " (A2+7y2A2)2(^A+C.)2 ' ^^•^> 

rM = = .(^Af (14) 

where p stands for the Cauchy principal value and J{lo) is the spectral density. Then, we have 



and 



TTJo (w - ?y A - + r(w)^ 



The integration in Eq. (|15p can be calculated numerically or approximately using residual theory. 



B. Measurement dynamics: quantum Zeno effect 

It is known that the quantum Zeno effect can effectively slow down the quantum decay rate of a quantum system. 
We study this effect using an approach that goes beyond the Markovian approximation and RWA. Here we consider 
the low-frequency bath as in Sec. II. A and derive the effective decay rate. The Hamiltonian H is given in Eq. (11), 

H w -i?7 Acr^ -f^WfcO^afc + ^77 A^^(a^cr_ +afccr+). (16) 

k k 

In this paper, the effective decay rate is defined in the same form as in Ref. [l^ . 
Write the wave function in the transformed Hamiltonian as 

i<i>(t))' = xWi;)iOfc)+E/3fcWit)iifc), (17) 

k 

with probability in the excited state at initial time |x(0)| • Substituting |$(t)) into the Schrodinger equation with 
H , we have 

.dx{t) _ A 



dt 2 



.dpk{t) ( 



^^xit) + J2Vkm, (18) 



z- 



dt 



uJk-^)Pk{t)+J2VkX{t)- (19) 



k 



When the transformations 



X{t)=x{i) exp 



77 A 
'i- — t 
2 



(20) 



Pk{t)^Pk{t) exp 



77 A 



are applied, Eqs. (IT51) and can be written as 

dm 



dt 
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(21) 
(22) 



dt 



= -iVkXit) exp [i(iLjfc - 77 A)t] 



Integrating Eq. and then substituting it into Eq. ([2^ . we obtain 



dm 

dt 



J x(t')exp[-z(wfc-?7 A)(t~t')] dt'. 



(23) 



(24) 



This integro-differential equation Eq. ([M)) is exactly soluble by a Laplace transformation. 

As for the present study of the quantum Zeno effect, i.e., using frequent measurements, it suffices to obtain the 
short-time behavior and the equation can be solved iteratively. With the initial excited-state probability amplitude 
x(0), in the first iteration, Eq. (|24|) is solved as 

m) 

k 



x(o) 



1 - /(* - E ^fe exp[-i(^fc - A)t'] 

■J 1. 



dt 



(25) 



We can approximately write x{t) using an exponential form: 



m 

x(o) 



exp — I {t — t) ^ exp —i{uJk — rj A)t 



dt 



exp < — < 



exp < —t 



, exp [—i {ujk — r] A)t] — 1 + i {ujk — rj A)t 



{^k - V 



1 2 sin 



,in {^^^^^t^ (^^. _ ^ A) t _ sin [(wfe - A) t] 



V 



t(uju~'n^y 



t{u}k-r] AY 



(26) 
(27) 

(28) 



The instantaneous ideal projections are assumed to be performed at intervals r. If single measurement, the probability 
amplitude is x{t = r). For a sufficiently large frequency of measurements, the survival population in the excited state 
is 



Poc(t = nr) = \xit = nr)p = |x(0)|^ exp[-7(T)i] 



(29) 



where the subscript "ee" refers to the initial and final excited state. And 7(t), with projection intervals r, is obtained 
as 



7(r) = 2^ du:Y.{^)\^^ 



2^kV 2sin^( 



k 

J{uj) 



uJk J T:{rjA — u))'^T ' 



— 2iT I duj 

Jo 4 



1 - 



uj — rj A 



1 2 



2sin"(2i^T) 
7r(77A — uj)^T 



(30) 



(31) 



cj + 77 A 

We now prepare the qubit to the dressed excited state exp[— S"] ||) |Ofe) (ct^ ||) = — |i)) at the initial time t — 0, which 
can be achieved by acting the operator on the ground state exp[— S*] |t) |Ofc) , 



|$(0)) = exph5] 1;) |Ofc) = a,exp[-5] |t) |Ofc) 



(32) 
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In this case, the initial state in the transformed Hamiltonian is \D \0k) and x(0) = 1- 

Note that in Eq. ((3T|) . the renormahzation factor ri of the characteristic energy A appears in the decay rate 7(t). 
This is different from the formulas for 7(r) in Refs. [2^ and (3]. In the case of spontaneous emission, the coupling 
strength between the electromagnetic field and atom is the fine structure constant 1/137, so it belongs to the weak- 
coupling case and rj then becomes extremely close to 1. Therefore, besides spontaneous emission, this result for the 
quantum Zeno effect can apply to other cases of strong coupling between the qubit and the bath. 

Under the RWA, the Hamiltonian becomes i?RWA as in Eq. (jG]). We prepare the initial excited state of H through 
the operator ax acting on the ground state under RWA |t) |Ofc) , <Tx |t) |Ofc) — \i) \0k)- Then, following the derivation 
in Refs. the decay rate is reduced to 

f°° V- /9fc\2 2sin2(A_±lr) , , 

,2/ A-tj_\ 



J(a;) 2sin^(^ 



= 27r / du; \' ' ^ (34) 

To compare with the high-frequency Ohmic bath, we choose the ordinary Ohmic bath with Drude cutoff: 

This is a realistic assumption for, e.g., an electromagnetic noise. In Eq. dSS]), is the coupling strength between 
the qubit and the Ohmic bath. The cutoff frequency Wc in the spectral density J°^(cj) is typically assumed to be the 
largest frequency in the problem. As for the low-frequency noise, we use JiiS) — 2au) j [uj"^ + A^), which is the same as 
in Eq. ([2|). The difference between these two baths is that A corresponds to an energy lower than the qubit energy. 
In the next section, we will show our numerical results for these two baths. 



III. RESULTS AND DISCUSSIONS 



To show the effects of either a low- or a high-frequency noise on the qubit states respectively, we study the dynamical 
quantities {(Jx{t)) and the quantum Zeno decay rate 7(t). The energy shift and the quantum Zeno decay rate exhibit 
evidently different features for the low- and high-frequency noises. Thus, these two quantities (energy shift and decay 
rate) can be used as criteria to distinguish the type of noise. In Fig. 1, we show the spectral densities J(a;) and J°^\lu) 
of the two baths in the cases of both weak dissipation a/A^ — 0.01 {a°^^ = 0.01) and strong dissipation otj ^ = 0.1 
{a°^ = 0.1). As examples, we choose A — 0.09A and A ~ 0.3A in Fig. 1(a) and 1(b), respectively. Here, each value 
of A corresponds to the position of the peak in the low-frequency spectral density. For an Ohmic bath, the cutoff 
frequency is fixed at uJc — lOA. As expected, these very different low- and high-frequency spectral densities should 
give rise to different decoherence behaviors of the qubit. In Fig. 1, we also show the characteristic energy of the 
isolated qubit A (see the vertical dotted line in Fig. 1). 



A. Non-measurement decoherence dynamics 



Results from the numerical integration of Eq. (26) are shown in Fig. 2. They are qualitatively consistent with the 
results obtained using residual theory. This indicates that the branch cuts considered in Refs. and do not affect 
the oscillation frequency. The time evolution of {(Jx{t)) is given in Fig. 2(a) for the case of weak coupling between the 
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FIG. 1: (Color online) The spectral density J(uj) of the low- and high-frequency baths, (a) The case of weak interaction between 
the bath and the qubit, where the parameters of the low-frequency Lorentzian-like spectrum are q/A^ — 0.01 and A — 0.09A 
(red solid curve), while for the high-frequency Ohmic bath with Drude cutoff the parameters are = 0.01 and luc ~ lOA 
(green dashed-dotted curve), (b) The case of strong interaction between the bath and the qubit, where the parameters of the 
low-frequency bath are q/A^ =0.1 and A = 0.3A (red solid curve) and the parameters of the high-frequency Ohmic bath are 
—0.1 and ujc = 10 A (green dashed-dotted curve). The characteristic energy of the isolated qubit is indicated by a vertical 
blue dotted line. Here, and in the following figures, the energies are shown in units of A. 



qubit and the bath, where a/A^ = 0.01 and A — 0.09A for the low-frequency noise and q;°'^ 0.01 and ujc — lOA for 
the Ohmic bath. Figure 2(b) presents the time evolution of {ax{t)) in the strong coupling case with the parameters 
a/A^ = 0.1 and A = 0.3A for the low-frequency noise, as well as = 0.1 and = lOA for the ohmic bath. As 
expected, the quantum oscillations of {(Jx{t)) dampen faster in the strong-coupling case. We approximately evaluate 
the oscillation frequency or the effective energy of the qubit, — ~ R{t^o) = 0, using the residue theorem. The 
decay rate can be obtained from r(Ld). We will now show the numerical values of rj in the corresponding cases. 
In Fig. 2(a), the renormalized factor rj = 0.98336, the oscillation frequency is ujq = 1.0225 A and the decay rate is 
r(a;o) 0.014654A for the low-frequency noise; while 77 = 0.98447, ojq = 0.97720A, and r(wo) 0.015318A for the 
Ohmic bath. In Fig. 2(b), the renormalized factor 77 = 0.91444, the oscillation frequency is ojq = 1.0868A and the 
decay rate is r(wo) = 0.11215A for the low-frequency noise, while t] = 0.84469, wq = 0.77221 A, and r(wo) = 0.13163A 
for the Ohmic bath. 

The energy spectral densities in Fig. 1 and the results in Fig. 2 indicate two opposite shifts of the characteristic 
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FIG. 2: (Color online) Time evolution of the coherence {(Txit)) versus the time t multiphed by the qubit energy spacing A. 
(a) The case of weak interaction between the bath and the qubit, where the parameters of the low-frequency Lorentzian-type 
spectrum are q/A^ — 0.01, A — 0.09A (red solid curve); while for the high-frequency Ohmic bath with Drude cutoff the 
parameters are — 0.01, ujc = lOA (green dashed-dotted curve), (b) The case of strong interaction between the bath and 
the qubit, where the parameters of the low-frequency bath are q/A^ — 0.1, A = 0.3 (red solid line), and for the high-frequency 
Ohmic bath are q"*" = 0.1, u)c = lOA (green dashed-dotted line). These results show that the decay rate for the low-frequency 
bath is shorter than for the high-frequency Ohmic bath. This means that the coherence time of the qubit in the low-frequency 
bath is longer than in the high-frequency noise case, demonstrating the powerful temporal memory of the low-frequency bath. 
Also, our results reflect the structure of the solution with branch cuts [l3|. The oscillation frequency for the low- frequency 
noise is ojo > A, in spite of the strength of the interaction. This can be referred to as a blue shift. However, in an Ohmic bath, 
the oscillation frequency is ujo < A, corresponding to a red shift. The shifting direction of the energy is independent of the 
interaction strength and only determined by the spectral properties. Thus, it can be used as a criterion for distinguishing the 
low- and high-frequency noises. 



energy A for the two kinds of baths considered here. These opposite energy shifts are equivalent to energy repulsion. 
The energy shift is determined by the interaction term. In the transformed Hamiltonian, the interaction is H[ in 
Eq. (|A.5p . Also, dipolar interactions, such as iJ/'-^ = g{a'^a- -f aa^) in the Jaynes-Cummings model, decrease the 
qubit's ground-state energy and increase its excited-state energy. Thus, now we ask the following questions: how a 
qubit is affected by either a multimode bath or a single-mode cavity? How a qubit is influenced by these two kinds 
of multimode baths: low-frequency and high-frequency ones? 
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For a low-frequency bath, the energy peak of the bath is located between the ground-state energy and the excited- 
state energy, that is to say the main part of the spectrum is in the region ujk < A. Then (as seen in Fig. 1) the 
interaction of tlic bath with the two qubit states is "opposite" , i.e. tlic ground-state energy becomes lower and the 
excited-state energy becomes higher. So the energy spacing for the case of a low-frequency bath exhibits a blue shift. 
This result is similar to the single-mode Jaynes-Cummings model. The energy difference of the two-state qubit is 
increased by the low-frequency bath. 

For a high-frequency cutoff Ohmic bath, the energy peak of the bath is located above the excited-state energy. 
So the main part of the spectrum is in the region Wk > A. The effect of the bath on the qubit mainly comes from 
the frequencies higher than the excited-state energy of the qubit. Then (as seen in Fig. 1) the bath repels both 
the excited-state and the ground-state energies to lower energies. But the effect of the high-frequency bath on the 
excited state is much larger than on the ground state. As a result, on the whole, the qubit energy difference in a 
high-frequency bath is red shifted. Thus, the effective energy difference of the qubit is reduced by the high-frequency 
bath. 

For example, if the initial state of the qubit is an excited state, in the interaction picture, the main part of the 
coupling is 

'^QkO-k exp(ia;A;t)cr_ exp(-iAi) 
fc 

= X]5fe«feC^- {cos[(wfe - A)t] +ism[{ujk - A)t]} , (36) 
fc 

where the real part contributes to the decay rate and the imaginary part results in the energy shift. For the low- 
frequency bath, the main part of the spectrum is in the region cok < A. Thus, the term for the energy shift is 

sin[(wfe - A)t] < 0. (37) 

However, for a high-frequency cutoff Ohmic bath, the main part of the spectrum is in the region oJk > A, where the 
term for energy shift becomes 

sin[(wfc - A)t] > 0. (38) 

These results show that the energy shift for the two kinds of baths moves in opposite directions. Note that there is a 
minus in the interaction term in the expressions for the dynamical quantities such as Eq. (21) and Eq. (33). The above 
observations help us understand why the energy levels repel. The contribution by the real part of the interaction on 
the decay rate will be discussed below, when studying the quantum Zeno effect. 



B. Measurement dynamics: quantum Zeno effect 

The quantum Zeno effect can be a useful tool to preserve the state coherence of a quantum system, with the help 
of repeated projective measurements. Below we investigate the quantum Zeno effect in the qubit system and propose 
another criterion for distinguishing low- and high-frequency noises. In general, without using the RWA, the effective 
decay rate can be obtained as 

.M.2,r^JHfl-^)'^f''^;'. (39) 
^' Jo \ uj + rjAj 7r(??A-w)2r ^ ^ 



11 



This expression includes three terms, i.e., the spectral density J{u>) of the bath, the projection time modulating 
function 

2sin2(IL^T) , , 

Tr{r] A — wfT 

and the interaction contribution function of both the rotating and counter-rotating terms 

The counter-rotating term contributing to /(w) is (w — r/ A)/ (w -|- A). If the RWA is applied, /(w) = 1. In Fig. 
3(a), the decay rate 7(t)/7o for the low-frequency bath is plotted in the weak-coupling case with a/A^ = 0.01. Here 
7o is the decay rate for |e) — ^ l^i) in the long-time limit with the RWA, 

7o=7(T-S^oo) = 27rJ(A)/4. (42) 

From the energy spectrum in Fig. 1(a), we can sec that 70 is proportional to the spectrum density of the bath, 
with the magnitiide corresponding to the crossing of the energy A of the isolated qubit and the bath spectrum. For 
comparison, we also plot 7rwa(t)/7o as a dashed-dottcd curve, with 



00 



J(u) 2sin^( 



2/ A- 



7RWA(r) = 2^/ duj \' ' ' ' . (43) 
Jo 4 7r(A - ujYt 

As wc know, 7(t)/7o < 1 means that repeated measurements slow down the decay rate 7(t) < 70, which is the 
quantum Zeno effect. On the contrary, 7(t)/7o > 1 means an anti-Zeno effect. The curves in Fig. 3(b) show results 
for the Ohmic bath with 0;°'^ = 0.01. In Fig. 4, we show the decay rate in the strong-coupling case for the low- frequency 
bath with a/A^ = 0.1 and for the Ohmic bath with 0!°*^ = 0.1. It can be seen that, for the low-frequency bath, the 
anti-Zeno effect appears as shown in Figs. 3 and 4. For the high-frequency cutoff Ohmic bath, the Zeno effect always 
dominates and no anti-Zeno effect occurs. Also we can see, from Figs. 3 and 4, that 7(t) and 7rwa(t) approach 70 
when the measurement interval r ^ 00. In particular, if r — >■ 0, F(w,r) — >■ 0. Thus, 7(r) 0. This implies that in 
a sufficiently short time interval of a projective measurement, the quantum Zeno effect occurs, regardless of the bath 
spectrum. When the interval r increases, the projection interval modulation function F(uj,t) displays a number of 
oscillations. Then, the energy peak of the bath spectriim will act on the decay rate and it is possible to implement the 
anti-Zeno effect. However, this result still depends on the function /(w) and the given spectral density J{iS) of the 
bath. Now, we emphasize again that the second term in the bracket of the function /(w) = [1 — (w — A)/ (w -|- A)]^ 
is due to the counter-rotating terms; when neglecting the counter-rotating terms, /(w) = 1. 

For low-frequency noise, the noise mainly comes from the region w < A, so /(w) > 1. Here, /(w) as well as F(w, r) 
magnify the effect of the energy peak of the bath spectrum. Thus, the counter-rotating terms accelerate the decay 
and the anti-Zeno effect occurs. 

In the high-frequency cutoff Ohmic bath, the noise mainly comes from the region a; > A, which leads to /(w) < 1. 
Thus, it is mainly the coimter- rotating term in /(cj) that reduces the effect of the energy peak of the bath on the decay. 
This slows down the decay and only the quantum Zeno effect can now take place. The projcction-intcrvals-modulating 
fmiction F{uj,t), together with the interaction-modulating function /(to), causes the Zeno effect to dominate in the 
high-frequency cutoff Ohmic bath. 
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a/A =0.01, A, = 0.09 
With RWA 




FIG. 3: (Color online) The effective decay 7(t)/7o, versus the time interval r between consecutive measurements, for a weak 
coupling between the qubit and the bath. In the horizontal axis, the time-interval r is multiplied by the qubit energy difference 
A. The curves in (a) correspond to the case of a low-frequency bath with parameters a/A^ = 0.01 and A = 0.09A (red solid 
curve), (b) corresponds to the case of an Ohmic bath with parameters = 0.01 and ujc ~ lOA (red solid curve). The green 
dashed-dotted curves are the results under the RWA when the same parameters are used. Note how different the RWA result 
is in (b), especially for any short measurement interval r. 



IV. SUMMARY 



In summary, we have studied a model of a qubit interacting with its environment, modeled either as a low- or 
as a high-frequency bath. For each type of bath, the quantum dynamics of the qubit without measurement and the 
quantum Zeno effect on it are shown for the cases of weak and strong couplings between the qubit and the environment. 
Our results show that, for a low-frequency bath, the qubit energy increases (blue shift) and the quantum anti-Zeno 
effect occurs. However, for a high-frequency cutoff Ohmic bath, the qubit energy decreases (red shift) and the 
quantum Zeno dominates. Moreover, for a high-frequency environment, the decay rate should be faster (without 
measurements) or slower (with frequent measurements, in the Zeno regime), compared to the low- frequency bath 
case. These very different behaviors of the quantum dynamics and the Zeno effect in different baths should be helpful 
to experimentally distinguish the type of noise affecting the qubit and protect the coherence of the qubit through 
modulating the dominant frequency of its environment. 
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FIG. 4: (Color online) The effective decay 7(t)/7o, versus the time interval r between successive measurements, for a strong 
coupling between the qubit and the bath. The time-interval r is multiplied by the qubit energy difference A. The curves in (a) 
correspond to the case of a low-frequency bath with parameters a/A^ =0.1 and A = 0.3A (red solid curve), (b) corresponds 
to the case of an Ohmic bath with parameters a"*" = 0.1 and ujc = lOA (red solid curve). The green dashed-dotted curves are 
the results under RWA when the same parameters are used. Note how different the RWA result is in (b), especially for any 
short measurement interval r. 
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Appendix: Examining the validity of the unitary transformation 



In this appendix, we show the main results of the canonical transformation to the qubit-bath Hamiltonian H 
considered here, and prove that the contribution of H2 to physical quantities is of the order Oi^g"^ and higher, so we 
ignore H2 in the calculations. We now apply a canonical transformation to the Hamiltonian H: 



F' = exp [S]Hc^Y> 



(A.l) 



with 



(A.2) 



Here a fc-dependent variable, = u)k/{ujk + 77 A), is introduced in the transformation. It is clear that this canonical 
transformation is unitary, because (exp [S\)^ = exp [~S]. The transformed Hamiltonian H' can now be decomposed 
in three parts: 



H' — Hq + H[ + H2, 



(A.3) 



with 



1 0^ 

k k 



(A.4) 



77 A - 



(A.5) 



where 



Hi, — — Adr < cosh 



A 



sinh 



ak] 



ry = exp 



E9k 
— 

k 



gfc 



(A.6) 



(A.7) 



Note that no approximation was used during the transformation, so H' = exp [S] iJ exp [—5*] is exact. Because the 
2 ___ 

constant term j^£,k{'^ — ^k) in Eq. (jA.4[) has no effect on the dynamical evolution, we neglect it. Considering 
at low temperatures, the multiple-step process is so weak that all the higher-order terms can be neglected. In the 
following derivation, we will prove that the contribution of H2 to physical quantities is of the order Ol^g"^ and higher, 
so we ignore H2 and obtain the effective transformed Hamiltonian H' = Hq + H[ . 

Now let us expand the first term of H2, cosh J2k Z^^kia^ ~ ^k) , as a series in gui^kl^k- We define Xk as Xk — 
gkik/^k, which is proportional to gk'- 



cosh 



9k_ 



/ ^ ^^k{al - ak) 



cosh 
1 f 



L k 

exp 



"^Xkia^ - ak 



exp 



Yxk{a+ - ak] 



(A.8) 
(A.9) 
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The first term in Eq. (IA.9[) can be written as 



exp 



t -ak) 



exp 



exp 



exp 



^Xka'l,^-Xkak 



(A.IO) 



where [^j, XkO-t ^ Sfc ^XfcO^fe] _ means the commutator of two operators. Using the commutation relation [a^, Ofe] — 
— 1, it is simplified to 



exp 



^Xkat,Y,-Xkak 



exp 



i \ " 2 

9 .Z^ 



exp 



.2 



E5fe J- 2 



(A.ll) 



which is the definition of 77. 

Afterwards, we expand Eq. (lA.lOp as follows, 



exp 



77 exp 



exp 



-Xkttk 



fe 
fe 



(Efe xfe( 



3! 



(-EfeXfeflfe)^ (EfeXfeOfc)^ 



2 3! 



(A.12) 
(A.13) 
(A.14) 
(A.15) 



l\, /c /c 

Now we see the second-order terms in Xk in Eq. (|A.15I) : 

^Xkaj^Xkak = ^Xfeflfeafe (diagonal) + ^ Xfea^Xfe'^fc' (off-diagonal). (A.16) 

fe fe fe fc^^fe' 

The off-diagonal terms are related to the multi-boson transition and their contributions to the physical quantities 
are fourth order in Xk- Furthermore, the initial state |Ofc) |t) is used in the calculation, so the direct effect of the 
second-order diagonal term on physical quantities is zero and its effect through the interaction will also be fourth 
order in Xk- Thus we now ignore the terms higher than second order in Xk in the following calculation and obtain 
the expression of Eq. (jA.10|) , 



exp 



^Xk{a+ - ttk) 



1 + E^'="fe ^^Xkttk 



Similarly, the second term in Eq. (|A.9p becomes 

E»(" 



exp 



fc - flfe) 



T] exp 



exp 



^Xkak 



^-^Xka^ +^Xkak 

k fe 



Therefore, Eq. (IA.9P is reduced to 

fffc 



cosh 



Y] —^k{a+ - Uk) 



— cosh 



L fe 
exp 



"^Xkiat -ak 
^Xk{a^ 



t-o-k) 



exp 



(A.17) 

(A.18) 
(A.19) 

(A.20) 

(A.21) 
(A.22) 
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In the same way, the third term in H2 is simphfied to 



sinh 



E 



fife 



6(4 



(A.23) 



In il2) W6 have subtracted the terms of zero and first-order in g^, which are included in Hq + H[, and only left the 
terms equal to and higher than second-order in gk, whose contribution to the physical quantities is of order 0[gfj 
and higher. Thus, H2 can be omitted. 

We expanded several series in the variable Xk, with the real variable rjAxk less than gk. That is: riAxk = 

rjAgk^.k/'jJk = V^9k/ (wfc + jyA) < gk- In other words, through this transformation we find a variable smaller than g^ 
for the series expansion. Therefore, our method can be extended to the case of strong interaction between the qubit 
and the environment. 
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